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Abstract The 16S rDNA amplicon high-throughput sequencing technique provides a robust and inexpensive approach 
to detect the gut microbiota of amphibians. Since different experimental protocols generate technical biases in drawing 
the gut microbiota profiles, the integrative analysis of gut microbiota produced by different studies must be performed 
with circumspection. In this study, we compared the efficacy of two DNA extraction methods (1.e., a phenol-chloroform 
method and TIANamp Stool DNA Kit) in describing intestinal and fecal bacterial communities of transplanted Asiatic 
toad (Bufo gargarizans) tadpoles. In terms of the DNA extraction quality De, DNA purity and yield rate) and the 
consistency in between fecal and intestinal microbiota structures (1.e., a and B diversity indices), the phenol-chloroform 


method was more robust than this commercial stool kit in profiling gut microbiota of tadpoles with feces. 
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1. Introduction 


Knowledge about the diversity and composition of 
gut microbiota has accumulated rapidly by culture- 
independent methods, especially the 16S rDNA amplicon 
high-throughput sequencing technique (Su et al., 2012). 
However, different experimental protocols concerned 
with material preparation, DNA extraction and PCR 
primer selection can generate technical biases in drawing 
the gut microbiota profiles (Larsen et al., 2015; Lozupone 
et al., 2013; Wagner Mackenzie et al., 2015). 

Since acquiring high quality genomic DNA is the 
prerequisite for the downstream analysis of gut microbiota 
structure, researchers have performed many studies on 
the effect of material types, sample storage conditions 
and DNA extraction methods (Choo et al., 2015; Ferrand 
et al., 2014; Larsen et al., 2015; Song et al., 2016; 
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Wagner Mackenzie et al., 2015). In spite of differences 
between intestinal and fecal microflora (Eckburg et al., 
2005; Nava et al., 2011), feces have been utilized as a 
noninvasive material to measure the intra- and inter- 
individual variation in composition and diversity of 
intestinal bacterial communities (Arumugam et al., 2011; 
Dethlefsen et al., 2008; Yatsunenko et al., 2012). Storage 
conditions for fecal samples could significantly influence 
microbiota profiles (Bahl et al., 2012; Choo et al., 2015; 
Maukonen et al., 2012), though some exceptions exist 
(Fouhy et al., 2015). DNA extraction protocols can 
produce technical variations due to different cell lysis and 
purification methods (Claassen et al., 2013; Maukonen 
et al., 2012; Wagner Mackenzie et al., 2015; Yuan et al., 
2012). However, these technical variations seem not 
to be large enough to distort the biological variations 
(Lozupone et al., 2013; Wagner Mackenzie et al., 2015). 
Recent research on gut microbiota has shed new light 
on the associations between gut microbiota and vertebrate 
physiology, development and evolution (Ley et al., 
2008; McFall-Ngai et al., 2013; Nicholson et al., 2012; 
Yatsunenko er al., 2012), while the number of intestinal 
microflora studies in amphibians remains far smaller than 
that in mammals (Bletz et al., 2016; Colombo et al., 2015; 
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Jiménez and Sommer, 2017; Kohl et al., 2013; Kohl et al., 
2014; Mashoof et al., 2013; Vences et al., 2016; Weng 
et al., 2016; Weng et al., 2017). Intensive studies of the 
gut microbiota in different taxa including amphibians are 
essential prerequisites to elucidate host-gut microbiota 
symbioses, e.g., phylosymbiosis (Brooks et al., 2016; Ley 
et al., 2008; Li et al., 2017a; Li et al., 2017b; Shapira, 
2016; Vences et al., 2016). In addition, the understanding 
of gut microbiota in amphibians could help us to take 
effective measures for amphibian conservation and 
cultivation (Jiménez and Sommer, 2017). Since many 
amphibians have been suffering from severe survival 
conditions (Hof et a/., 2011), noninvasive approaches will 
always be an optimal choice for the study of symbioses 
between amphibians and gut microbiota. However, it 
remains unclear whether the feces can be applied as 
a noninvasive material to study the gut microbiota of 
amphibians. 

In this study, we aimed to test the efficacy of a phenol- 
chloroform method and a commercial fecal reagent kit 
(TIANGEN Biotech Co., Ltd.) in describing intestinal 
bacterial communities of Asiatic toad (B. gargarizans) 
tadpoles by feces. Specifically, DNA extraction quality of 
different methods and sample types was tested in terms 
of three parameters, 1.e., A260/A280, A260/A230 and 
DNA yield rate. Furthermore, the structural consistency 
between bacterial communities was evaluated based on 
16S rDNA amplicon high-throughput sequencing. 


2. Materials and Methods 


2.1. Transplantation of B. gargarizans eggs In 
February-March 2016, we sampled seven broods of B. 
gargarizans eggs from Xinyang City of Henan Province 
in China (Table S1). A ~10 cm-length chalaza was taken 
from each brood of eggs, and then was hatched in labs 
using plastic cylinders (1 L) filled with about 0.5 L 
dechlorinated drinking water. Subsequently, about 15 
larvae in each brood were reared together with boiled 
green vegetable leaves rich in cellulose or fish foods 
rich in protein. We applied semi-natural conditions for 
rearing all larvae, i.e., a water change and food feeding 
per three days without controlling the light, humidity 
and temperature. All procedures used in this study were 
approved by the Animal Care and Use Committee of 
Xinyang Normal University. 


2.2. Preparation of intestinal and stool samples When 
tadpoles developed into the lower limb stage (Gosner 
35-40), we collected intestinal and stool samples. We first 
collected 10 samples (S1—S10) of mixed fecal sediments 


from cultivation water by using aseptic injectors (Table 
S1). These mixed stool samples spin-dried in a centrifuge 
were applied for comparing DNA extraction quality of 
the phenol-chloroform method and the TIANamp Stool 
DNA Kit. In addition to the DNA extraction quality, 23 
tadpoles (S11—S33) were further selected to compare 
the microbiota structure between the phenol-chloroform 
method and the commercial stool kit (Table S1). Tadpoles 
(S11-S16 and S22—S27) were individually cultivated for 
less than 24 hours after a water change. Then we got their 
spin-dried stool samples through the same approach for 
mixed stool samples. We sacrificed tadpoles using 7596 
ethanol solution before the extraction of the gut samples 
into sterile microcentrifuge tubes. All intestinal and stool 
samples were stored in a —20?C freezer before the DNA 
extraction. 


2.3. Protocols of DNA extraction and DNA quality 
evaluation Phenol-chloroform method.—The 500 uL 
SDS (196) and each weighed sample were mixed in a 
sterile 1.5 mL microtube and bathed in the water for 10 
min at 60?C. During the water bath, microtubes were 
overturned and blended for three times. Subsequently, 
30 uL EDTA (0.5 M) and 20 uL protease K (20 mg/mL) 
were added into each microtube, and the mixture was 
bathed in the water for 1 hour at the same temperature of 
the previous step. And then the samples were centrifuged 
at 12 000 rpm for 5 min at 4?C. The supernatant 
fraction was transferred to a new 1.5 mL microtube, and 
blended with an equal volume of phenol-chloroform- 
isoamyl alcohol (25:24:1 in volume). The samples were 
centrifuged at 12000 rpm for 5 min at 4?C once again. 
The supernatant fraction was transferred to a new 1.5 
mL microtube, and blended with an equal volume of 
chloroform-isoamyl alcohol (24:1 in volume). The 
samples were centrifuged at 12000 rpm for 5 min at 4?C 
once again. The supernatant fraction was transferred 
to a new 1.5 mL microtube, and blended with twofold 
absolute ethanol. After cooled at —20?C for 20 min, the 
samples were centrifuged at 12 000 rpm for 10 min at 
4°C. After the removal of the supernatant fractions, the 
DNA precipitates were washed in 200 uL ethanol (70%) 
for 5 min, and then the mixtures were centrifuged at 12 
000 rpm for 5 min at 4?C. This step was repeated three 
times. The dry DNA was dissolved in 40 uL TE buffer 
solution. 

TIANamp stool DNA kit.—TIANamp Stool DNA 
Kit simplifies the DNA isolation by a fast spincolumn 
procedure. The protocol recommended by the 
manufacturer was utilized to extract the metagenomic 
DNA. 
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DNA quality evaluation.—The A260/A280, A260/ 
A230 and DNA yield rate were applied for the DNA 
quality evaluation. The A260/A280, A260/A230 and 
concentration of DNA products were determined by using 
NanoVue Plus Spectrophotometer (GE Healthcare Inc., 
Germany). The DNA yield rate was given by the ratio of 
the DNA concentration to the sample weight. The DNA 
products extracted from feces or intestines of tadpoles 
(S11-S16 and S22-S27) were subsequently stored in a 
—20?C freezer until the high-throughput sequencing of 
16S rDNA amplicons. 


2.4. 16S rDNA amplicon high-throughput sequencing 
The library construction of 16S rDNA amplicons and 
high-throughput sequencing on MiSeq (Illumina Inc., 
USA) were achieved in a commercial company (Genergy 
Inc., China). Specifically, the hypervariable regions of 
V3-V4 in the bacterial 16S rDNA were amplified from 
the microbiota DNA products using the universal primer 
pair 341F-CCTACGGGNGGCWGCAG and 785R- 
GACTACHVGGGTATCTAATCC (Klindworth et al., 
2013). The amplicons were generated by a two-step, 
tailed PCR on the DNA products in terms of the 16S 
Metagenomic Sequencing Library Preparation protocol 
with some modifications (http://support.illumina.com/ 
downloads/16s metagenomic sequencing library ` 
preparation.html). The volume of each PCR solution 
was 25 uL, which consisted of 10 ng DNA, 2.5 uL 10 x 
Takara Ex Taq Buffer, dNTP (2.5 mmol/L each), 5 pmol/ 
L forward primer, 5 pmol/L reverse primer, 0.1 uL Takara 
Ex Taq, and ddH,O. The 1* PCR condition was 94°C for 
3 min, 20 cycles (94?C for 10 s, 55?C for 15 s, 72?C for 
30 s) and 72°C for 7 min. The 2" PCR condition was 
94°C for 3 min, 5 cycles (94°C for 10 s, 55°C for 15 s, 
72°C for 30 s) and 72°C for 7 min. Finally, the 16S rDNA 
amplicons quantified by Quant-iT™ PicoGreen® dsDNA 
Assay Kit (Life Technologies, USA) were paired-end 
sequenced on the Illumina MiSeq platform. 


2.5. Bioinformatic analyses of 16S rDNA amplicons 
The raw paired-end reads, after the removal of barcodes, 
were filtered by Trimmomatic software with three 
minimum thresholds (Bolger et al., 2014), i.e., terminal 
base quality score equal to 25, average quality score 
equal to 25 in sliding windows of 50 bp with a step of 
1 bp, fragment length equal to 100 bp. The paired-end 
sequences were merged with the flash software for a 
minimum overlap of 10 bp and a maximum mismatch 
proportion of 0.2. The merged sequences including 
ambiguous bases were excluded. Subsequently, the 16S 
rDNA reference database downloaded from NCBI was 


Profiling Gut Microbiota of Amphibian Tadpoles by Feces 3 


utilized for the validation of V3—V4 regions. In this study, 
16S rDNA amplicons were successfully sequenced in 32 
samples excluding S14 S, S17 G and S20 G (Tables 1, 
S1). 


Table 1 Sequencing results for successfully sequenced 32 samples. 


Parameter Value 
Number of total sequences 811311 
Number of total bases 3.74 x 108 
Minimum sequence length 124 
Maximum sequence length 510 
Mean of sequence lengths 460.71 
Median of sequence lengths 465 
Standard deviation 11.81 
GC percentage 0.5384 
N50 465 


Operational taxonomic units (OTUs) were assigned 
to the 16S rDNA amplicons using UPARSE software 
with the identity threshold of 97% (Edgar, 2013). The 
annotation of representative OTU sequences was executed 
through sequence alignments on the RDP database using 
the RDP classifier (confidence threshold = 0.8) (Wang 
et al., 2007). 

The mothur program was used to calculate a diversity 
indexes (Schloss et al., 2009), i.e., richness (OTU 
observed, ACE and Chaol), diversity and evenness 
(Shannon, Simpson and Shannoneven), and Good’s 
coverage. The rarefaction and Shannon curves were 
drawn to measure whether the sequencing depth was 
appropriate for the richness and diversity calculation. To 
statistically analyze the differences in a and f diversity 
indexes between samples, 16 881 Ge, minimum number 
of sequences in all samples) sequences in each sample 
were subsampled. In addition, we drew stackbars and 
Venn charts to show the taxonomic composition and 
abundance of samples. Furthermore, we used Past 
software (v3.14) to execute principal coordinates analysis 
(PCoA) based on Bray-Curtis dissimilarity calculated 
from taxonomic abundances to analyze the f diversity 
of microbiota structure (Hammer er al., 2001). The 
transformation exponent was set to the default value De, 
2). 


2.6. Statistical analyses of the DNA extraction quality 
and the microbiota structure To compare the DNA 
extraction quality of phenol-chloroform method and the 
TIANamp Stool DNA Kit, we utilized SPSS software 
(v20.0, IBM Corporation) to perform Mann-Whitney 


4 Asian Herpetological Research 


Vol. 9 


Rank Sum Test and two-way ANOVA (two factors: 
experimental material and DNA extraction method) in 
terms of A260/A280, A260/A230 and DNA yield rate. 
Pairwise multiple comparisons were performed in two- 
way ANOVA with Bonferroni correction when overall 
significance level was less than 0.05. In addition, fitting 
curves, which sketched the relationship of the DNA yield 
rate and sample weight, were compared between the 
phenol-chloroform method and the commercial stool kit. 
To test the effect of experimental material and DNA 
extraction method on the microbiota structures, we 
executed two-way ANOVA to compare a diversity 
indexes of intestinal and fecal microbiota in S11—S33 
tadpoles. We hypothesized that the factors, i.e., host 
genetic background and diet, could bias the comparative 
analysis of microbiota structure. Therefore, we used 
paired ¢-Test or Wilcoxon Signed Rank Test to compare 
a diversity indexes of gut microbiota between 10 pairs of 
littermate tadpoles (S11—S16, S18—S21 and S24-833), 
i.e., SII. G-SI6 G, S18 G-S21 G versus S24 G-S33 G 
(Gut ph vs Gut kit). Subsequently, the comparisons of a 
diversity indexes between intestinal and fecal microbiota 
were performed on 12 tadpoles (S11-S16 and 822-827), 
i.e., SII S-S16 S versus S11 G-S16 G (Stool ph vs 
Gut ph) and S22 S-S27 S versus S22 G-827 G (Stool 
kit vs Gut. kit). The intestinal/fecal metagenomic DNA of 
S11-S16 and S22-S27 was extracted using the phenol- 
chloroform method and TIANamp Stool DNA Kit, 
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Kit Phenol-choroform 


A260/A280 


respectively. 

In addition to permutational multivariate ANOVA 
(PERMANOVA), we also performed a Mantel test on the 
Bray-Curtis dissimilarities calculated from taxonomic 
(De, OTU, genus and phylum) abundances in the above 
three cases, 1.e., Gut ph vs Gut kit, Stool ph vs Gut ph 
and Stool kit vs Gut kit. The number of permutations 
was set to 9999 for PERMANOVA and Mantel tests. To 
test which taxonomies significantly affected the structural 
divergence of microbiota in these three cases, Lda 
Effective Size (LEfSe) were executed with recommended 
options (http://huttenhower.sph.harvard.edu/galaxy/). 


3. Results 


3.1. Efficacy comparison in terms of DNA purity and 
yield rate Among the tests of the between-method 
difference in A260/A280, A260/A230 and DNA yield 
rate, a significant difference was detected in A260/A230 
(Mann-Whitney Rank Sum Test: U = 205, n, = 23, n, = 
27, P = 0.04) and DNA yield rate (U = 66, n, = 15, n, 

17, P = 0.02) but not A260/A280 (U = 300.5, n, = 23, n, 
= 27, P = 0.85) (Table S1; Figures 1, 2). Whereas when 
the two-way ANOVA was applied (Table 2), a significant 
between-method difference was detected in A260/A230 
and DNA yield rate but not in A260/A280. A significant 
between-material difference was detected in A260/A280 
and A260/A230 but not in DNA yield rate. No significant 


Kit Phenol-choroform 


A260/A230 


Figure 1 Comparisons of A260/A280 and A260/A230 between TIANamp Stool DNA Kit (n = 23) and phenol-chloroform method (n = 27). 


"P «0.05. 
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Table 2 Two-way ANOVA on A260/A280, A260/A230 and DNA yield rate. 


Source of variation Digiesot Sum of squares Mean square F P 
freedom 
Materials 2 0.14 0.07 3.73 0.03 
A260/A280 Methods 1 0.06 0.06 3.09 0.09 
Materials x Methods 2 0.07 0.03 1.83 0:17 
Materials 2 1.94 0.97 7.22 < 0.01 
A260/A230 Methods 1 0.70 0.70 5.21 0.03 
Materials x Methods 2 0.63 0.31 2.33 0.11 
Materials 2 48 055.72 24 027.86 1.85 0.18 
DNA yield rate Methods 1 123 332.77 123 332.77 9.49 < 0.01 
Materials x Methods 1 5490.65 5490.65 0.42 0.52 


interaction was detected between these two factors (i.e., 
experimental material and DNA extraction method). In 
the pairwise multiple comparisons, A260/A230 (avg. 
+ std. err.) of stool samples (1.05 + 0.11, n = 12) was 
significantly different from that of gut samples (1.51 4 
0.08, n = 23) and mixed stool samples (1.51 + 0.10, n = 
15). Furthermore, the fitting curves for DNA yield rate 
and sample weight show that the DNA yield rate in the 
phenol-chloroform method and the commercial stool kit 


decreases along with saturated DNA extraction materials 
in solutions (Figure S1). The phenol-chloroform method 
has a greater DNA yield rate than the commercial stool 
kit does at a given sample weight. 
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3.2 Efficacy comparisons in terms of microbiota 
structures In Gut ph vs Gut _ kit, total number of OTUs 
was 121, among which 4 belonged to Gut ph (3.3%), 4 
belonged to Gut Kit (3.3%) and 113 belonged to between 
groups (93.4%) (Figure 3). In Stool ph vs Gut ph, total 
number of OTUs was 100, among which 6 belonged 
to Stool ph (6%), 8 belonged to Gut ph (8%) and 86 
belonged to between groups (86%). In Stool kit vs Gut | 
kit, total number of OTUs was 114, among which 9 
belonged to Stool kit (7.9%), 7 belonged to Gut kit 
(6.1%) and 98 belonged to between groups (86%). The 
Venn chart indicates that the majority of OTUs could be 
extracted by using the phenol-chloroform method or the 


Phenol-choroform 


Figure 2 Comparison of DNA yield rate between TIANamp Stool DNA Kit (n = 17) and phenol-chloroform method (n = 15). 'P < 0.05. 
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commercial stool kit from intestines or stools. 

No significant differences in o diversity indexes were 
detected in between- material groups and between-method 
groups by using two-way ANOVA (Table S2). When 
host genetic background and diet factors were taken into 
account, we detected no significant differences in almost 
all o diversity indices in Gut ph vs Gut kit, Stool ph vs 
Gut ph and Stool kit and Gut kit (Table 3; Figure 4). One 
exception was that the richness values (i.e., number of 
observed OTUs) showed a significant difference between 
Stool kit and Gut kit (P — 0.02, n — 6) (Table 3; Figure 
4). 

None of the Bray-Curtis dissimilarities calculated 


from OTU, genus and phylum abundances showed a 
significant between-group difference in Gut ph vs Gut ` 
kit (Table 4; Figures 5, S2, S3). However, only the Bray- 
Curtis dissimilarity based on genus abundances gave 
rise to a weak between-group correlation (R — 0.32, P 
= 0.04). In the PERMANOVA for Stool ph vs Gut ph 
(Table 4; Figures 5, S2, S3), the OTU, genus and phylum 
abundances also showed an insignificant between-group 
variation in terms of Bray-Curtis dissimilarity. In addition, 
both Bray-Curtis dissimilarities calculated from OTU and 
genus abundances showed a significant between-group 
correlation, i.e., R = 0.64 (P = 0.04) and R = 0.58 (P = 
0.02). As for the Stool kit vs Gut kit (Table 4; Figures 5, 


Table 3 Alpha diversity indexes (avg. + std. dev.) calculated from OTU tables rarefied to 16 881 sequences per sample in three cases, i.e., 
Gut ph vs Gut. kit (n = 9), Stool ph vs Gut ph (n = 5) and Stool kit vs Gut kit (n = 6). 


Gut ph vs Gut. kit 


Stool ph vs Gut ph Stool kit vs Gut kit 


Gut ph Gut. kit Stool ph Gut ph Stool kit Gut. kit 
Richness 62.33 + 11.56 60.67 + 9.45 63.60 + 4.72 58.20 + 12.28 69.17 + 8.86 58.17 + 12.83 
Shannon 1.84 + 0.63 2.08 + 0.55 1.82 + 0.52 1.62 + 0.52 1.76 + 0.49 1.80 + 0.46 
Chaol 75.00 + 12.73 70.14 + 10.95 71.40 + 5.49 69.02 + 11.62 78.4 £11.48 — 68.035 15.12 
ACE 78.52 + 14.58 71.90 + 11.26 73.02 + 3.91 72.65 + 16.75 79.81+11.95 66.78 + 12.90 
Simpson 0.33 + 0.20 0.25 + 0.16 0.32 + 0.18 0.39 + 0.18 0.37 + 0.19 0.30 + 0.17 
Coverage 1.00 + 0.00 1.00 + 0.00 1.00 + 0.00 1.00 + 0.00 1.00 + 0.00 1.00 + 0.00 
Shannoneven 0.44 + 0.14 0.50 + 0.12 0.44 + 0.12 0.40 + 0.11 0.41 € 0.11 0.44 + 0.10 
Gut ph Stool kit 
x 
€ 


Figure 3 The distribution of OTUs in fecal and intestinal microbiota based on two DNA extraction methods. 
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Table 4 PERMANOVA and Mantel test on the Bray-Curtis dissimilarities calculated from taxonomic (Oe, OTU, genus and phylum) 
abundances in three cases, i.e., Gut ph vs Gut. kit (n = 9), Stool ph vs Gut ph (n = 5) and Stool kit vs Gut. kit (n = 6). 


Gut ph vs Gut. kit 


Stool ph vs Gut ph Stool kit vs Gut. kit 


OTU genus phylum OTU genus phylum OTU genus phylum 

PERMANOVA  Permutation No. 9 999 9 999 9 999 9 999 9 999 9 999 9 999 9999 9999 

F 1.25 1.02 0.50 0.63 0.44 0.58 4.33 4.37 2.51 

P 0.22 0.40 0.67 0.66 0.80 0.54 0.02 0.02 0.12 

Mantel test Permutation No. 9 999 9 999 9 999 9 999 9 999 9 999 9 999 9999 9999 

Correlation (R) 0.16 0.32 -0.06 0.64 0.58 0.28 0.65 0.64 —0.17 

P 0.21 0.04 0.54 0.04 0.02 0.22 0.01 0.001 | 0.75 
100 S2, S3), both Bray-Curtis dissimilarities calculated from 
90 OTU and genus abundances also showed a significant 
80 between-group correlation, i.e., R — 0.65 (P — 0.01) and 
70 R — 0.64 (P — 0.01). Nevertheless, significant between- 
60 group variations was detected in both Bray-Curtis 
50 "un dissimilarities based on OTU and genus abundances 
40 BGut_kit — (PERMANOVA: F,,,—4.33, P = 0.02) and F; 1) = 437, P 

30 = 0.03). 

20 In Gut ph vs Gut kit, we identified by the LEfSe 
10 analysis that three genera of phylum Proteobacteria 
M — -=r showed significant between-group divergence, i.e., 
90 Sphingorhabdus biased to Gut_ph, and Coxiella 
80 biased to Gut_kit along with an unnamed genus in 
70 order Rhizobiales (Figure 6). As for Stool ph vs Gut 
60 H Stool ph ph, five genera (i.e., Hydrogenophaga, Rhizobium, 
50 Gut. ph Brevundimonas, an unnamed genus in order Rhizobiales, 


0 
100 


Stool kit 


Gut kit 


Figure 4 Comparisons of a diversity indices in three cases, i.e., 
Gut ph vs Gut kit, Stool ph vs Gut ph and Stool kit vs Gut _ 
kit. The bars and error bars represent mean values and standard 
deviations, respectively. 'P « 0.05. 


an unnamed genus in family Sphingomonadaceae) in 
phylum Proteobacteria biased to Stool ph, and genus 
Clostridium XlVa in phylum Firmicutes biased to 
Gut ph (Figure 6). In addition, significant between- 
group divergences of four phyla (i.e., Proteobacteria, 
Bacteroidetes, Verrucomicrobia, Firmicutes) composed 
of 15 genera were detected in Stool kit vs Gut kit, 1.e., 
12 genera in the phyla of Proteobacteria, Bacteroidetes 
and Verrucomicrobia biased to Stool kit, and 3 genera 
in the phyla of Proteobacteria and Firmicutes biased to 
Gut kit (Figure S4). The shared between-material bias 
of genera in Stool ph vs Gut ph and Stool kit vs Gut kit 
included four genera in phylum Proteobacteria biased to 
Stool (i.e., Hydrogenophaga, Rhizobium, Brevundimonas, 
an unnamed genus in family Sphingomonadaceae) and 
one genus in phylum Firmicutes biased to Gut De, 
Clostridium XIVa). 


4. Discussion 


The phenol-chloroform method is a classic and cost- 
effective approach to extract eukaryotic or prokaryotic 
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Figure 5 Principal coordinates analysis (PCoA) based on Bray-Curtis dissimilarity calculated from OTU abundances in three cases, i.e., 
Gut ph vs Gut kit, Stool ph vs Gut ph and Stool kit vs Gut kit.. 
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Figure 6 LEfSe results in the case of Gut_ph vs Gut_kit and Stool_ph vs Gut_ph. 


genomic DNA (Ausubel, 2002). The DNA yield rate of 
the phenol-chloroform method is significantly higher than 
those of several commercial methods (Wagner Mackenzie 
et al., 2015). Similarly, this study demonstrated the DNA 
extraction efficiency of the phenol-chloroform method 
outweighed that of TIANamp Stool DNA Kit in terms 
of DNA yield rate. In addition, the phenol-chloroform 
method produced a better A260/A230 than TIANamp 
Stool DNA Kit did, which probably resulted from 
more effective DNA washing in the last step of phenol- 
chloroform method. Due to the poor operationality and 
technical overlook of mucosal bacteria in separating gut 
contents from intestines, we extracted the mixture of gut 
microbiota DNA and host genomic DNA from intestinal 
samples. Nevertheless, we applied the universal primer 
pairs to specifically amplify V3—V4 regions of bacterial 
16S rDNA (Klindworth er al., 2013). It has been reported 
that the mechanical treatments (e.g., bead beating) of 
samples can produce more efficient cell lysis of Gram- 
positive bacteria and yield high amounts of DNA (Ferrand 
et al., 2014; Guo and Zhang, 2013). However, no 
method has become a gold standard for 16S rDNA high- 
throughput sequencing (Yamagishi et al., 2016). Methods 
with a mechanical treatment face a tradeoff between cell 


lysis efficiency and DNA disruption level. Due to the easy 
digestion of tadpole intestines and feces we did not take 
the mechanical treatments into account. 

The composition and diversity of gut microbiota are 
undoubtedly affected by multiple factors, e.g., host genetic 
background, dietary profile and environmental situation 
(Dabrowska and Witkiewicz, 2016; Davenport, 2016; Jin 
et al., 2017; Voreades et al., 2014). The annotation on 
the gut microbiota structure can be susceptibly biased by 
the experimental protocols (Choo et al., 2015; Claassen 
et al., 2013; Yuan et al., 2012). However, many studies 
have demonstrated that biological variations outweigh 
technical variations generated by DNA extraction 
methods (Salonen er al., 2010; Wagner Mackenzie 
et al., 2015; Wesolowska-Andersen er al., 2014) and 
sample storage conditions (Blekhman et al., 2016; Fouhy 
et al., 2015). Here we homogenized and minimized 
biological variations in between groups as far as possible. 
Nevertheless, it is impossible to eliminate each interfering 
factor. For instance, the littermate tadpole pairs used in 
Gut ph vs Gut kit possess a similar but not an identical 
genetic background, thereby the between-method 
heterogeneity is possibly enhanced. From the analyses 
of Gut ph vs Gut kit, the technical variation generated 
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by DNA extraction methods was outweighed by the 
inter-subject variation. However, the phenol-chloroform 
method and the commercial stool kit probably resulted in 
a significant inconsistency in the structural composition 
of microbiota, e.g., OTU and phylum abundances. To 
ensure biological differences outweigh systematic biases, 
we had better use the identical standardized protocols for 
the comparative analysis of gut microbial consortia. 

Although the microbiota structure shows a significant 
variation between feces and intestinal contents in the case 
of the sophisticated intestines (Gu er al., 2013), feces 
has been applied as an effective noninvasive material for 
the study of gut microbiota in mammals. Larsen et al. 
(2015) used ribosomal intergenic spacer analysis to reveal 
that the microbiota structure of fishes was similar but 
significant different in feces and intestines. However, their 
preparation procedure of fecal and intestinal samples, 
i.e., the fecal samples were squeezed from the intestinal 
samples, probably enhance the between-material 
variations. In this study, the fecal microbiota based on 
the phenol-chloroform method can more efficiently 
reflect the gut microbiota in terms of composition and 
diversity. On the contrary, the feces and intestines possess 
more inconsistent microbiota structures deduced from 
the commercial stool kit. When we applied the PCoA to 
compare the samples fed with different food types, the 
experimental material factor seemed to be dominant rather 
than genetic background and food type factors in the case 
of commercial stool kit. We will explore and discuss the 
effects of multi-factors (e.g., genetic background and 
food type) on the microbiota structure of Asiatic toad 
tadpoles in a further study. Even though no significant 
difference between fecal and intestinal microbiota was 
detected in the phenol-chloroform method, between- 
material variations do exist, e.g., four genera in phylum 
Proteobacteria were more abundant in feces and one 
genus in phylum Firmicutes in intestines. We argue that 
the inconsistency possibly resulted from the moderate 
interference of microorganisms in water and tadpole skin 
to feces. 


5. Conclusion 


According to the DNA extraction quality and structural 
comparisons between fecal and intestinal microbiota, 
the phenol-chloroform method is probably more robust 
than a commercial fecal reagent kit (TIANGEN Biotech 
Co., Ltd.) in evaluating the gut microbiota structure 
of amphibian tadpoles with feces. To the best of our 
knowledge, this study provides the first evidence that 


feces of amphibian tadpoles can be applied as an effective 
noninvasive material for the study of gut microbiota. 
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Table S2 Two-way ANOVA on a diversity indexes. 


Source of variation df F P 
Richness Materials 1 0.92 0.34 
Methods 1 0.83 0.37 
Materials x Methods 1 1.86 0.18 
Shannon Materials 1 0.00 0.97 
Methods 1 0.07 0.79 
Materials x Methods 1 0.02 0.88 
Chaol Materials 1 1.44 0.24 
Methods 1 0.84 0.37 
Materials x Methods 1 3.85 0.06 
ACE Materials 1 1.53 0.23 
Methods 1 0.36 0.56 
Materials x Methods 1 4.10 0.05 
Simpson Materials 1 0.09 0.77 
Methods 1 0.05 0.82 
Materials x Methods 1 0.20 0.66 
Coverage Materials 1 1.93 0.18 
Methods 1 0.03 0.86 
Materials x Methods 1 3.23 0.08 
Shannoneven Materials 1 0.05 0.82 
Methods 1 0.02 0.90 
Materials x Methods 1 0.16 0.69 
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Figure S1 Relationship between DNA yield rate and sample weight in TIANamp Stool DNA Kit (n = 17) and phenol-chloroform method (n 
= 15). 
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Figure S2 PCoA on Bray-Curtis dissimilarity calculated from genus abundances in three cases, i.e., Gut_phvs Gut_kit, Stool ph vs Gut ph 
and Stool_kit vs Gut_kit. 
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Figure S3 PCoA on Bray-Curtis dissimilarity calculated from phylum abundances in three cases, i.e.,Gut_ph vs Gut. kit, Stool ph vs Gut ph 
and Stool kit vs Gut Kit. 


PS Gutkit — Stool kit Stool kit vs Gut kit 


U 1 Li 
Bacteria Proteobacteria. Betaproteobacteria. Burkholderiales 
Bacteria Proteobacteria. Betaproteobacteria. Burkholderiales.Comamonadaceae 


Bacteria. Proteobacteria. Betaproteobacteria. Burkholderiales. .Comamonadaceae:Hydrogenophaga 


Bacteria. Proteobacteria. Betaproteobacteria. Burkholderiales. Comamonadaceae. Hydrogenophaga. Other 

` Bacteria. Proteobacteria. Gainmaproteobacteria.Alteromonadales 

: Bacteria] Proteobacteria {Gammaproteobacteria Alteromonadales.Shewanellaceae 

: Bacteria.| Proteobacteria. Gammaprotéobacteria. Alteromónadales. Shewanelláceae. Shewanella 

Bacteria. Proteobacteria.Gammaproteobacteria Alteromonadales.Shewanellaceae/Shewanella. Shewanella. | Xiamenensis 
: Bacteria. Bacteroidetes. Bacteroidia. Bacteroidalés. Bacteroidaceae 


Bacteria, Bacteroidetes.Bacteroidia.Bacteroidales.Bacteroidaceae.Bacteroides 
Bacteria. Bacterdidetes. Flavobacterita. Flavobacteriales. Flavobacteriaceae 

: Bacteria. Bacteroidetes. Flavobacteriia 
Bacteria.Bacțeroidetes.Flavobacțeriia „Flavobacterialps.Flavobacteriaceąe-Flavobacterium 

: : Bacteria. Bacteroidetes. Flavobacteriià. Flavobacteriales 

‘Bacteria. Bacteroidętes. Flavobacteriia Flavobacteriales.| Flayobacteriaceae. Flavobacterium. Other 
` Bacteria. Proteobacteria Alphaproteobacteria: Rhi zobiales Rhizobiaceae: Rhizobium.Other 


Bacteria. Proteobacteria.Gammaproteobacteria.Alteromonadales.Shewanellaceae.Shewanella.Other 
Bacteria. Proteobacteria Betaprotdobacteria. Burkholdériales. Comamonadaceae. Acidovorax 


Bacteria Proteobactéria Betaproteobacteria | Burkholderiales:Comamonadaceae Acidovorax. Other 
Bacte ia.Proteobacteria. Alphaproteobacteri .Rhizobiales.| Rhizoblaceae. Rhizobium 
Bacteria Protéobacteria.Alphaproteobacteria Rhizobiales Rhizobiaceae 


Bacteria Proteobacteria.Betaprotepbacteria.Burkholderiales. Burkholderiales i incertae sedis 

Bacteria, Proteobacterid Alphaproteobacteria. .Caulobacterales.t Caulobacteraceae 

Bacteria Proteobacteria, Alphaproteobacteria. Caulobacterales. Caulobacteraceae.Brevundimonas 

H Bacteria. Proteobacteria/Alphaproteobacterià. Caulobacterales 

Bacteria.Proteobacteri iphaproteobacteria aulobacterales.Caulobacteraceae.Brevundimonas.Other 


Bacteria.Proteobacteria. .Alphaproteobacteria.Sphingomonadales 


: Bacteria.Proteobacteria.Alphaproteobacteria.Sptiingomonadales.Sp! lingomonadaceae 
Bacteria.Proteobacteria ;Betaproteobacteria;Burkholderiales.Comamonadaceae.Limnohabitans.Other 
f Bacteria. Proteobacteria Alphaproteobacteria. Rhodobacterales. Rhodobacteraceae. Other 
Bacteria. Proteobacteria. Alphaproteobacteria. Rhodobacterales. Rhodobacteraceae. Other.Other 


Bacteria. roteobacteria.| Betaproteobacteria. Burkholderiales. Burkhdlderiales ; incertae . Sedis.Sphaerotilus 
Bacteria.| Proteobacteria. Betaproteobacteria. Burkholderiales. Burkholderiajes ; incertae sedi Sphaerotilus. Other 
Bacteria. Bakteroidetes. Cytophagia. Cytophagales 

;Bacteroidetes.Cytophagia.Cytophagales.Cytophagaceae 

` S Bacteria Bacteroidetes.Cytophagia 
Bacteria.Proteobacteri “Alphaproteobactéria. Sphingomonadsles.Sphingomonadeiceae.Other. Other 


Bacteria Proteobacteria Alpha protepbacteria -Sphingomionadales.Sphingomjonadaceae.Other 


Batteria. Bacteroidetes. Cytophagia.! Cytoptiagales. Cytophagacéae. Flectobacillus 
Bacteria. Bacteroidetes. Cytophagia. Cytophagales. Cytophagaceae. Flectobacillus: uncultured | Bactergidetes t bacterium 


Bacteria. Verrucomicrobia Vefrucomicrobiae. Verrucomicrobiales. Vertucomicrobiaceae. Akkermansia.Other 
Bacteria. Proteobacteria. etaproteobacteria, Burkholderiales.t Oxalobacteraceae. Undibacterium.Other 
Bácteria. Proteobactetia. Betaproteobactetia. Burkholderiales.Oxalobacteraceae 


Bacteria.Proteobacteria | Betaproteobacteria| Burkholderiales. Oxalobacteraceae.Undibacterium 


Bacteria. Firmicutés. Clostridia.Clostridiales. ER Eed 


ER 
EJ 
v 
RU 
o 
PLE 
D 
9 
a 
E 
a 
E 
D 
I 
a 
a 
US 
[3 
LÀ 
ki 
Ed 
8 
9 
k3 
S 
O 
e. 
D 
o 
[e] 
9 
B 
EI 
Fn 
am 
4 
B 
E 
e$. 
e 
& 
s 
2 
S 
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Figure S4 LEfSe results in the case of Stool kit vs Gut Kit. 


